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Abstract 

We study the Landau free energy for a uniaxial ordering, taking into account 
two Umklapp terms of comparable strengths (those of the third and fourth 
order). Exploring the analogy with the well-known nonintegrable classical 
mechanical problem of two mixed nonlinear resonances, we complete the pre- 
vious studies of the corresponding phase portrait by calculating numerically 
periodic solutions, including those far from the separatrices. It is shown that 
in the physical range of parameters only periodic configurations are absolutely 
stable. We determine for the first time the complete thermodynamic phase 
diagram and show that, in contrast to some earlier claims, the wave number 
of the ordering does not pass through the devil's staircase, but through a fi- 
nite number of steps which decreases as the amplitudes of the Umklapp terms 
increase. 
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The sine-Gordon model [T|,0] which is widely used in the studies of the uniaxial 
incommensurate-commensurate (IC) transitions ||, is based on a few intricate approxi- 
mations. By one of them only one, presumably dominant, phase dependent Umklapp term 
is retained in the Landau expansion of the free energy. In this Letter we critically examine 
this approximation by including two Umklapp terms which favor different commensurabil- 
ities, i.e. those of the third and fourth order. Like in other extensions of the sine-Gordon 
model |§^§, one encounters a nonintegrable problem, well-known in the modern classical 
mechanics PJlOU . As a consequence, the thermodynamic phase diagram for the IC transition 
qualitatively differs from that of the simple sine-Gordon model. 

It is well-known that continuous nonintegrable models |7],|| and discrete models |4]-|6| 
for the IC transition lead to similar descriptions of uniaxial orderings. The analyses were 
however usually made for limited ranges of values of physical parameters which enter into 
models. Also, the continuous models were examined only in the parts of phase space close 
to separatrices, characterized by dilute soliton lattices. With such restrictions it was not 
possible to determine complete thermodynamic phase diagrams which are the basis for any 
discussion of particular real systems. In this Letter we combine two numerical methods, with 
the aim to investigate the entire phase space and to find the regions of absolute stability for 
particular solutions in the plane spanned by two control parameters (out of three present 
in the model). It will be shown that, due to the gradual disappearance of quasiperiodic 
solutions and the formation of chaotic layers, the phase diagram is (almost) entirely covered 
by periodic configurations. 

The model is defined by the free energy functional 

F= [ dx[-A 2 + B cos (30 + 3(- -5Ax) +C cos (40-4^)1 (1) 
J 2 ax 6 

with cf) being the phase of the order parameter. 5^ = Q — Q4, where Q is the wave number 
of ordering in the absence of Umklapp terms and Q4 = 27r/4 is the commensurate wave 
number of the fourth order. (We put everywhere the lattice constant equal to unity.) By 
assumption Q is placed somewhere between Q 4 and the commensurate wave number of the 



third order Q3 = 27r/3. In that case the Umklapp terms of third and fourth order which 
are included into the Landau expansion ([[]) might be of comparable strengths. E. g. the 
Umklapp term of m-th order for charge and spin density waves is proportional to the m-th 
power of the amplitude of the order parameter Thus, at temperatures not far below 
the critical temperature T c of the incommensurate ordering, the third order Umklapp term 
should not be omitted due to the dominance of the coupling constant B with respect to C. 
On the other hand, for small values of 64 (e. g. for £4 <C tt/6 — 64) the closeness of the wave 
number to the commensurate point Q4 favors the fourth order Umklapp term. The similar 
arguments can be invoked for other types of uniaxial orderings. Note that the example ([[]) 
is the most interesting one since it includes two Umklapp terms which are of the lowest 
possible orders in the amplitude of the order parameter. (The model for the commensurate 



lock-in at Q = tt is of essentially different type [|T2|JT^ 

The competition between two periodicities in the expansion ([I]) takes place in the non- 
linear terms which are formed on the purely symmetry grounds. It is present in both the 
dependence and the " nonautonomous" x-dependence of these terms. We recall that the 
starting functional of Fradkin et al |7]], which is appropriate for the approximate treatment 
of the ordering in the presence of mobile defects contains only the competition of the 
latter type. 

The extremalization of the functional (H) leads to the Euler-Lagrange (EL) equation 

0" + 3B sin [30 + 3 (tt/6 - 8 4 )x] + AC sin (40 - A5 A x) = 0. (2) 

Among all solutions of this equation we are looking for that which is the absolute minimum 
of the mean free energy density in the eq. ([!]) To this end we develop a procedure 



which is complementary to those usually performed in the classical mechanical analyses of 
the phase portrait of eq. (fj) (see e. g. refs. [§,|K]]). Standard methods of direct integration 
of eq. (0), which were widely used in the calculations of quasiperiodic and chaotic orbits, 
cannot be extended to the isolated periodic solutions, which are orbitally unstable in the 
most interesting cases of phase portraits with finite chaotic layers. However, it appears that 



just these solutions have the main role in the thermodynamic phase diagram. We therefore 
calculate them separately, by treating eq. (|2]) as the boundary value problem. Combining 
these two methods we can follow all relevant trajectories in the phase portrait and extract 
the trajectory for which the functional ([I]) has the absolute minimum. 

The functional (|I|) and the EL equation (0) do not possess the continuous translational 
symmetry in x, but are invariant under the combinations of discrete translations in both x 
and fl5l| . All solutions of eq. (H) can be thus grouped into degenerate subsets of solutions 



connected by these translations. In the further considerations of the periodic solutions we 
enumerate these subsets and represent each subset by one solution. 

It is easy to determine possible periods of the periodic solutions, defined by 

<j)(x + P) = <j>{x) + (j) P . (3) 

Inserting eq. (|3|) into the EL equation (0), one gets the necessary conditions 

P = 4k + 31, <f)p = 5 4 P - Ztt/2, (4) 

where k and I are integers. The values of P are <pp are thus discretized (in contrast to the 
simple sine-Gordon case) and denote allowed intervals of variables x and <j) i n the boundary 
value problem. Note that the solutions with finite values of <pp are rotational {i.e. have 
finite winding numbers). 

The boundary conditions for the periodic solutions with the values of P and <pp given 
by eq. (|J) are defined by <fiki( x = Xi) , 4>ki( x = Xf = Xi + P) = 4>ki( x — x i) + 4>p an d 
4>' k i(x = Xi) = 4>' k i(x = Xf). Using two different numerical methods |L6| we undertake a 
systematic search of such conditions for periods P up to ~ 100. Both numerical procedures 
give at most two solutions with the same values of k and /, defined by two isolated points 
in the two-dimensional plane of boundary conditions [(fiki(x = x i),(p'ki( x = x i)\ 1FH - We 



noticed |T5[] that these points sometimes have jumps and discontinuities as the functions 
of the parameters B and C. In particular, a given periodic solution may cease to exist 
as the nonlinearity strengthens, the general tendency being that the solutions with smaller 
periodicities are more robust as B and/or C increase. 



When both exist, two solutions with same values of P and <pp have different values of 
the free energy ([[]). Both have odd symmetry [the only possible for rotational solutions], but 
with different positions of inversion points [xq,4>(xq)]. It easily follows from the differential 
equation (fj) that xq necesarily has only integer and half-integer values, while 4>{xq) = 
Xo^4 — mir/A, m = 0, ±1, ±2, .... This is indeed the case with all our numerical solutions. 
With the odd symmetry taken into account, the numerical procedure is further simplified, 
since it remains to determine only one continuous boundary condition [4>' k i(x = Xo) = 4>' k i(x = 
xq + P)} and two discrete numbers [xq and ^{xq)] per each solution. This analysis |15| shows 



that the particular values of xq and 4>(xo) for a given solution are related to the evenness 
and oddness of the corresponding integers k and /. Further below we follow only periodic 
solutions with the lower free energy ([!]), and characterize each such solution by a pair of 
integers (k,l). 

The structure factor for a given periodic solution (k, I) has equidistant satellites, whose 
positions with respect to the the wave number Q are given by 

<f) P 2mr . 2n-l/2 n , , , n 

^ = T + ^ = 5i + 7T ^kTk> n = 0,±l,±2,... (5) 

The position of the fundamental satellite go is determined by the total slope of the function 
<f>ki(x), while the higher harmonics q n ,n ^ appear due to its periodic modulation. 

In the next step we determine mean free energies of particular solutions. For illustration, 
we show in Figs.la-c mean free energies of periodic, quasi-periodic and chaotic solutions 
with one fixed = 3/2) = 0] and one variable [<p' Q = 4>{xi = 3/2)] initial condition, 

choosing 5^ = it/ 12 and three representative values of B and C. For all periodic solutions 
shown in these figures Xi = 3/2 coincides with the point of odd symmetry Xq. Since these 
figures remain qualitatively same for other choices of the "initial" point Xj, they enable some 
general conclusions on the relationship between free energies of various types of solutions. 

The situation in which the parts of the phase space with the rotational flow are still 
free of chaos is presented in Fig. la. It is evident that there is a periodic solution which 
has lower mean free energy than all other calculated periodic and quasiperiodic solutions. 



The same numerical conclusion is obtained for other finite values of B and C. It should be 
however mentioned that the comparison of mean free energies of periodic and quasiperiodic 
configurations becomes more and more subtle as B and/or C decrease, so that in this limit 
one cannot rigorously exclude the absolute thermodynamic stability of some quasiperiodic 
solution 

We turn from this intricate and numerically unpleasant limit to the intermediate values 
of B and/or C for which both, quasiperiodic and chaotic, layers are present in the rotational 
part of the phase space [Fig. lb], and finally to the regime || 

4B + VC > tt/12 (6) 

in which there is no quasiperiodic solution and this part is covered by a single chaotic 
layer (Fig.lc). [Note that the wings on the left and right sides in Figs.l.a-c come from 
the thermodynamically uninteresting oscillating solutions of eq. (Q).] Since the numerical 
integration of eq. (|2|) cannot extract a single chaotic solution (due to its orbital instability), 
the plateaux in Figs. lb, c strictly represent mean free energies within a given chaotic layer, 
i.e. each plateau represents one isolated chaotic layer in the phase space. Still, it is plausibly 
expected |7| that each particular chaotic solution fully diffuses inside "its" layer, so that 
its mean free energy coincides with one of the plateaux in Fig.l. Periodic solutions are 
immersed into these chaotic layers as isolated points. Their free energies are discontinuous 
(needle-like) minima well below all chaotic plateaux [see Figs. lb (inset) and lc]. As it was 
already noticed, the number of surviving periodic solutions decreases with the increase of 
"chaoticity" of the phase space, measured by the left-hand side of the inequality @. 

The numerical results illustrated by Figs.la-c are the basis for the construction of the 
thermodynamic phase diagram in the three-dimensional space of parameters B, C and S4. 
One possible visualization is presented in Figs. 2 and 3 in which we fix B, vary C and 
#4 |18|], and follow the fundamental wave number go [eq. @] and the free energy of the 
configuration with the lowest value of the functional (|I|). For C large enough [e. g. such 
that the inequality @ is satisfied], the system passes through small number of periodic 
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configurations as 64 varies. Other existing periodic solutions have free energies which are 
larger and increase as the period P increases ||15|| . Since it does not seem probable that 
the further solutions, which have even higher periods and were not calculated numerically, 
would not follow this tendency, we conclude that there are no other solutions which would 
fill the phase diagram between the neighboring plateaux shown in the Fig.2. 

As C decreases the number of periodic configurations increases through a series of bi- 
furcations, while the dependences of the free energy and the wave number on £4, although 
discontinuous, approach that of the simple sine-Gordon limit (7=0. The number of ab- 
solutely stable periodic configurations still remains finite (and smaller than the number of 
numerically calculated periodic solutions), even for C ~ 0.001 — 0.002 when there are no 
chaotic layers in the rotational part of phase space (except those in the regions of separa- 
trices) (Fig. la). In other words, our numerical results suggest that the staircase from Fig.2 



remains "harmless" |TJ] {i.e. not fractal) as long as C is finite (or at least not too small). 
The harmless staircase with first order transition between neighboring periodic solutions 
is expected to be closely connected with the participation of dense soliton lattices in our 
phase diagram. Then it is inappropriate to use the picture of well separated solitons with a 
repulsive effective interaction, which was usually the basis for the conclusions that the wave 
number go may pass through the complete or incomplete devil's staircase 

In conclusion, we stress that the periodic solutions are of the central importance for the 
thermodynamics of systems modelled by the free energy ([!]). In the physical range of pa- 
rameters (B and C finite), only these solutions (i.e. the commensurate configurations) have 
the property of absolute thermodynamic stability. In the phase space they are isolated and 
immersed in a mostly chaotic environment, whose mean free energies represent characteris- 
tic energetic barriers which the system has to overcome in order to replace one periodicity 
by another. Such barriers might be intrinsic sources ]20| of memory effects and thermal 
hystereses which are frequently observed in uniaxial systems with a variable wave number of 
ordering [Q]. We recall that the thermodynamic significance of periodic and chaotic solutions 
is just opposite to that of the classical mechanical counterpart of the model (|l|). Mechanical 



periodic trajectories cannot be realized due to their orbital instability, so that the rotational 
motion of a strongly nonintegrable classical mechanical system is entirely chaotic. 
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FIGURES 

FIG. 1. Mean free energy vs <fa for S 4 = tt/12, B = 0.002, C = 0.002 [a], B = 0.008, C = 0.006 
[b] and B = 0.02, C = 0.02 [c] of periodic (A), quasiperiodic (O) and chaotic (solidO) solutions. 
The (k, I) indices for the periodic solution with the lowest free energy (solidA) are (3, 4) in the 
figures [a] and [b], and (1,1) in the figure [c]. The insets in figures [a] and [b] are enlarged 
neighborhoods of the free energy minima. 

FIG. 2. The fundamental wave number qo vs C and 84 for B = 0.02. The fractions at 
particular steps denote the ratios 2(^4 — (ft))/ 71 " = ?/(4fe + 3£). In front of the dotted cross-section 
the rotational part of the phase space does not contain quasiperiodic solutions but only a single 
chaotic layer and isolated periodic solutions [see eq. (§)]. 

FIG. 3. Free energy of the absolutely stable periodic configurations vs C and 5^ for B = 0.02. 
The fractions and the dotted cross-section have the same meaning as in Fig. 2. Note the reversed 
scales at the 64 axes in Figs. 2 and 3. 
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